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Abstract. 


3y  use  of  the  multinle-time-scale  method,  the  low  density  expansion 
is  carried  to  the  order  of  the  triple  collision  integral.  The  validity 
of  3ogoliubov's  assumption  that  the  multiple  distribution  depends  function¬ 
ally  on  a  single  particle  distribution  is  carefully  examined.  It  is  found 
that  such  an  assumption  is  valid  except  locally  for  those  particles 
which  have  a  large  senaration  at  a  time  t  and  which  have  their  relative 
velocity  so  oriented  that  they  were  in  collision  at  t  =  0,  Since  this 
local  breakdown  is  very  selective,  the  trinle  collision  integral  which 
is  found  in  the  literature  is  still  correct.  As  a  by-oroduct  of  the 
multinle-time-scale  method,  the  rate  at  which  a  system  approaches  the 
kinetic  state  is  obtained;  it  is  also  found  that  up  to  the  order  we  have 
considered  the  Maxwellian  distribution  is  the  only  solution  at  thermal 


equilibrium. 


I.  Introduction. 


During  the  past  decade,  various  methods  have  been  developed  for  the 
derivation  of  the  Boltzmann  collision  integral  from  first  principles  of 
mechanics  (Bogoliubov,1  Prigogine  and  his  co-vorkers ,2  Green, ^  and 
Kirkwood  and  his  co-workers*).  Recently  the  correction  to  the  Boltzmann 
(binary)  collision  integral  has  also  been  obtained  by  various  investi¬ 
gators  (Choh  and  Uhlenbeck,^  M,  Green, ^  Rice,  Kirkwood  and  Harris,"^ 

Q 

and  Resibois  ).  The  relationship  of  various  methods  was  first  discussed 
9  10 

by  Cohen  *  who  established  within  certain  assumptions  the  equivalence 

of  the  method  of  Bogoliubov  and  Green.  More  recently  Sandri11  and 

12 

McCune,  Sandri  and  Frieman  studied  the  problem  of  the  low  density  ex¬ 
pansion  of  the  B-B-G-K-Y  hierarchy  by  the  method  of  multiple-time-scales. 
They  found  that  the  triple  collision  integral  obtained  by  the  nrevious 
investigators  is  divergent.  The  physical  nature  of  this  divergence  was 
not  studied  in  any  detail.  On  the  other  hand,  the  nature  of  the  higher 
order  corrections  to  the  Fokker-Planck  collision  integral  was  carried 
out  in  great  detail  by  Su.1^  The  nature  of  local  singularities  in  the 
solution  for  correlation  function  was  displaced  there.  Well-behaved 
correction  terms  to  the  Fokker-Planck  collision  integral  were  obtained. 

In  this  paper,  we  shall  use  the  method  of  multiple-time-scales  to 
investigate  the  nature  of  Boltzmann  collision  integral  and  its  correction 
terms  under  the  assumption  of  spatial  homogeneity.  This  systematic  ex¬ 
tension  procedure  offers  a  simple  way  to  demonstrate  to  what  extent  the 
Bogoliubov  assumption  regarding  the  higher  distribution  functions  as 
functionals  of  the  first  distribution  is  valid.  It  is  found,  up  to  the 
triple  collision  level,  in  contrast  to  the  works  in  Refs,  11  and  1?,  that 


Bogoliubov's  assumption  is  valid  except  for  a  special  local  region  in 

two-particle  phase  space,  i.e.,  for  two  particles  having  a  large  separation 

at  time  t  and  with  their  relative  velocity  oriented  in  such  a  way  (say 

in  reference  to  their  relative  position  vector)  that  they  were  in  collision 

at  t  =  0.  This  same  local  breakdown  was  first  found  in  the  investigation 

of  a  weak  coupling  gas"^  which  leads  to  Fokker-Planck  equation  in  its 

lowest  order  approximation. 
lU 

Cohen,  in  his  investigation  of  the  triple  collision,  also  found 
it  necessary  to  neglect  the  kind  of  singularity  mentioned  above  in  order 
that  the  equation  for  the  first  order  two-particle  function  would  admit 
any  solution  at  all.  He  remarked  that,  for  the  "large  majority"  of 
cases,  the  equation  for  the  two-particle  function  due  to  the  triple 
interaction  was  solvable  and  has  a  well-behaved  solution.  He  seemed 
also  to  encounter  the  same  difficulty  in  Ref.  10  where  he  suggested  a 
coarse-graining  procedure. 

We  shall  indeed  see  that  the  secularitv  we  have  found  can  be  elimi¬ 
nated  by  an  averaging  process  in  the  momentum  space,  provided  we  restrict 
ourselves  to  the  triple  collision  level.  As  one  goes  to  higher  orders 
such  an  averaging  is  no  longer  sufficient  to  smooth  out  the  local  singu¬ 
larity. 

It  is  suggested  by  the  nature  of  the  singularity,  that  a  new  scaling 
be  applied  to  the  hierarchy  equation  in  this  special  region  in  phase 
space.  The  unit  of  time  and  length  in  this  region  is  taken  to  be  the 
mean  free  time  and  mean  free  path  respectively  in  contrast  to  the  duration 
of  collision  and  range  of  interaction  outside  the  special  region.  There 
is  no  longer  a  separation  of  time  scales  in  this  region.  The  hierarchy 


equations  have  to  he  solved  Jointly  on  one  time  scale  (kinetic  time 
scale).  However,  we  note  that  the  one-particle  function  is  affected 
by  the  two-particle  function  only  "grossly,"  i.e.,  through  the  following 
integral  operator: 


Here  x  is  the  relative  position  between  particles  1  and  2,  and  <(>  is  the 
interparticle  potential.  The  subscript  1  (or  2)  refers  to  the  particle 
1  or  2.  The  velocity  integration,  above  orovides  the  coarse-graining 
process  which  we  mentioned  earlier;  this  is,  however,  not  sufficient  at 
higher  orders  because  the  singularity  in  the  two-particle  function  then 
becomes  too  strong.  Since  in  the  special  region,  the  separation  of  the 
two  particles  has  to  be  large,  if  the  range  of  the  pair  potential  is  very 
limited,  the  contribution  of  any  singularity  at  large  x  will  be  strongly 
de-emphasized  by  the  potential  force  di»/dx  under  the  integral  sign.  We 
therefore  conclude  that  even  though  Bogoliubov's  assumption  becomes 
locally  invalid  in  a  local  region,  the  triple  collision  integral  obtained 
in  the  literature  is  correct.  The  procedure  of  solving  the  hierarchy 
equations,  such  as  that  given  by  Choh  and  Uhlenbeck,"’  can  apparently  be 
pushed  to  higher  order,  provided  the  pair  potential  between  the  particles 
is  sufficiently  short-ranged. 

We  shall  formulate  our  problem  as  an  initial  value  problem.  Since 
the  relevent  correlations  between  particles  are  those  created  through 
the  interaction  of  particles,  the  initial  values  for  all  correlation 
functions  will  be  taken  to  be  zero.  The  effects  of  the  initial  values  of 
correlation  functions  on  the  one-particle  kinetic  equation  were  studied 


in  Fef.  ij.  Within  the  framework  of  the  multi nle-time-scale  formulation, 
all  the  physical  chanres  in  the  problem  anoenr  in  their  own  expropriate 
time  scales  automatically:  relaxation  towards  the  kinetic  stage  on  the 
fast  time  scale  (of  the  order  of  the  collision  time),  kinetic  evolution 
on  the  slow  time  scale  (of  the  order  of  the  mean  free  time).  For  the 
spatially  homogeneous  system  which  we  shall  analyze,  the  time  scales 
end  automatically  at  the  kinetic  time  scale,  which,  of  course,  is  what 
one  would  exnect  physically , 


II.  Binary  Collision  Integral. 

The  first  three  members  of  the  B-B-G-K-Y  hierarchy  for  a  spatially 
homogeneous  system  under  the  low  density  approximation  are  as  follows: 


3f 

at 


dx_  dv 


d£  .  3£_ 

■2  dx  3v 


(1) 


If  ♦  V.„  .  it  .  i|t  .  (-f-  -  T-)  IrU)  ♦  f(l)f(2)  ] 
at  —12  a jc  m  ax  av^  av^  — 


(2) 


n 

e  — 

m 


d£  di3  H  •  tgT"  [f(l)e<l-*.23)  ♦  h(*,i,123)] 


-  [f(2)g(x-C,13)  ♦  h(x,x-£,123)]> 


3h  3h  _  „  3h 

—  ♦  v  •  — ►  ♦  V  • 

at  -12  ax  -13  n 


1  3b,  .  ,  j>_  _a_ 

n  ac.  *  av,  '  a v3 


x  ao(x-c) 

c  h 
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-3 

leading  term  of  oraer  e. 


Here  f  is  the  one-particle  distribution  function;  g  and  h  are  the  two- 
and  thrce-nart ic le  correlation  functions  which  are  relnted  to  the  tvo- 
and  three-narticle  distribution  functions  as  follows: 

f(l)  =  F1(l) 

g( 1?)  =  r2(2)  -  f(l)f(P) 

h( 123 )  =  F3(123)  -  f{ l)g(12)  -  f(2)*(l3)  -  f(3)g(l2) 

The  time  variable  in  the  above  functions  is  omitted  for  simplicity.  The 
numbers  within  each  bracket  indicate  the  snecific  particles  in  which  ve 
are  interested. 

x_  and  are  the  position  vectors  of  narticle  1  relative  to  particles 
2  and  3  respectively.  The  small  parameter  e  equals  the  average  number  of 
particles  within  the  range  of  the  pair  potential  which  is  assumed  to  be 
renulsive  and  finite  in  range. 

The  time  and  length  in  Eqs,  (l)  to  (3)  have  been  normalized  bv  the 
duration  of  a  collision  and  the  range  of  the  pair  potential  respectively. 
The  mean  free  path  of  close  collisions,  in  which  we  are  interested,  can 
be  given  as  follows: 


where  n  is  the  average  density  and  r^  is  the  radius  of  the  collision  cross 
section.  Identifying  the  latter  as  the  range  of  the  oair  potential,  we 
see  that  the  ratio  of  the  range  of  notential  to  the  mean  free  oath  (or 
the  ratio  of  the  collision  time  to  the  mean  free  time)  is  e.  ror  very 
small  e,  we  have  therefore  two  distinct  time  scales  in  our  problem:  a 
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fast  time  which  is  of  the  order  of  the  collision  time  and  a  slow  time 

which  is  of  the  order  of  the  mean  free  time.  In  the  method  of  nultiole- 

time-scale,  the  slow  time,  as  well  as  the  fast  time,  is  treated  as  in- 

* 

dependent  variables,  i.e.,  we  extend  f(t)  to  f(t,  Et,  ezt,  ...).  The 
extra  freedom  introduced  is  used  to  demand  that  the  solution  for  f  be 
well-behaved  in  the  limit  of  t  ,  i.e.,  no  secular  behavior  for 
g(t  -»“>).  Formally  this  extension  in  the  time  variable  is  equivalent  to 
the  expansion  of  the  time  derivative  in  the  following  way: 


3_  . 
3t 

where 

dt 

_ 0 

dt 


3  ^  3 

3t  £  3t  + 
o  1 


and  so  on. 


(10 


We  now  expand  f,  g,  and  h  in  simple  power  series  in  e,  for  example, 


f  =  f 


(o) 


+  ef 


(1) 


♦  e2f 


(2) 


(5) 


From  the  zero  order  equations  of  f,  it  is  seen  that  f^°^  is  t^  in¬ 
dependent.  The  equation  for  (x, 12,tQ )  is 


3t 

o 


♦ 


—12 


3x  m  3x 


3  '  (o)  _  1 
3v^  *  m  3x 


1|-)f(o)(l)f(o)(2) 

(6) 


Its  solution  is  easilv  found  (since  f^°^  is  t  independent)  as 

o 


rio> (x,12,t  )  *  [3  .  (12)  -  li  f(o)(l)f(o)(2) 
o 


(7) 


The  inclusion  of  tOt  and  so  forth  is  purely  a  mathematical  procedure. 
A  discussion  of  the  termination  of  tire  scales  is  given  at  the  end  of 
t:-.e  next  section. 


where  (12)  =  exp  [-tQ  H( 12 )  ]  is  the  two-particle  streaming  operator  and 
o 

H  ( 12?)  -  v  .  1.  _  i  .  (J _ L.) 

-12  3x  m  3x  l3v  3v2' 


The  one  particle  function  up  to  this  order  of  approximation  is 


3f ( 1 )  .  3f^o)  n 


3t 


at. 


(o) 


d*. . 


dx  dv  -  - 

—  —2  dx  3v 


(8) 


If  we  integrate  this  equation  in  t  ,  the  second  term  on  the  left  side 

o 

of  Eq.  (8)  will  be  proportional  to  tQ.  To  avoid  the  secularitv  due  to 
3f^°Vat  in  the  solution  for  f^  on  the  time  scale  t^,  we  require 

^  d-2  S  *  37"  1S-»(12)  (9) 


3f 


(o) 


at. 


17 

m 


and 

(l2>  -  =_U2>Ji-(0)li)f(0)U>  (10) 

o  —  — 1  o 

The  integral  term  in  (9)  can  be  transformed  into  the  usual  Boltzmann 
collision  integral  (binary).  It  is  seen  from  Eq.  (9)  that  the  kinetic 
evolution  is  in  a  time  scale  of  order  1/e  compared  with  the  duration  of 
a  collision.  It  was  shown  that  the  former  is  .lust  the  mean  time  between 
collisions . 

Equation  (10)  gives  the  transient  towards  the  kinetic  stare  governed 

bv  Eq.  £-;).  The  (t  )  time  behavior  of  the  source  function  on  the  right 
•  o 

side  of  Eq.  (10)  can  be  analyzed  as  follows.  The  pair  potential  has  teen 
assumed  to  have  a  finite  range.  At  time  to,  because  of  the  factor 
d«/dx  in  Eq.  (10),  the  particle  2  must  be  within  the  range  of  particle 


1  in  order  that  the  integral  be  different  from  zero,  ^or  sufficiently  large 
t  ,  S  (12)  will  definitely  bring  particle  2  outside  the  range  of  particle 

O  •  v 

*  0 

1.  If  the  pair  potential  has  an  exnonential-like  tail,  the  source 

function  of  Eo.  (10)  will  decay  exnonentiallv  in  t  as  t  tends  to  in- 

c  o 

finitv.  Therefore,  in  that  case,  the  approach  to  the  kinetic  stage  is 
exponentially  fast,  in  the  time  scale  of  the  collision  time. 


III.  Triple  Collision  Integral. . 

The  zero  order  three-narticle  correlation  function  can  be  obtained 
from  ( 3 )  as 


h(o)(x,^,123,t  )  »{[S  (123)  -  1] 

1  O 

o 


-  [S_t  (12)  -  1] 

o 

-  [S_t  (23)  -  1) 

o 


(11) 


-  [S_t  (31)  -  1])  f(o)(l)f(o)(2)f(o)(3) 
o 

where  S  (123)  is  the  three-particle  streaming  operator,  exp  [-tQ  H(l23)]t 
o 

and  ii(l23)  is  the  three-particle  Hamilton  operator: 


H(l?3) 


a  a  i  a* 

V  •  — ♦  V  •  — ^ 

-1?  ax  —13  r.  ax 


a  r  ^1  a  4> 

3Vp  m  3^ 


1  a  b(x-£) 
B  “IT 


# 


The  pair  potential  had  been  assumed  to  be  repulsive. 
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The  integrand  of  this  integral  in  the  limit  of  t  •*  »  is  nlotted 

o 

in  Pig.  1.  The  integral  on  the  second  line  of  (13)  is  given  by  the 
shaded  area  in  the  praoh.  If  x  is  finite,  the  collision  of  particles 

n 

1  and  2  (if  they  hannen  to  collide)  will  occur  at  a  finite  f  for  any 
nonvanishing  v_jp,  The  above  mentioned  area  is  therefore  finite  and 
hence  is  the  integral.  However  if  we  let 


x 


<T  — ip'o 


(lb) 


the  integral  becomes  singly  infinite,  where  v_  are  constants  with 
a  >  1  and  ^  related  to  the  range  of  the  pair  potential.  En.  ( iU )  in¬ 
dicates  a  very  careful  selection  of  the  points  in  the  nhase  space  of  the 
two  interacting*  particles.  We  want  the  narticles  to  have  a  lar^e 

separation  at  time  t  and  to  have  their  relative  velocity  oriented  in 
o 

such  a  way  (in  reference  to  the  direction  of  their  relative  nosition 

vector,  sav),  that  they  were  in  collision  at  t  =0.  T'uch  secularitv  was 

o 

also  found  in  the  weak  courling  expansion. ^  Mathematically  i t  is 
clear  that  the  secularity  is  due  to  the  infinite  range  of  1?)  in 

the  special  region  of  phase  space  mentioned  above.  We  observe  that  the 
solution  for  g^°^(12,to)  given  by  t'n.  (7)  cannot  be  valid  in  this  special 
region,  since  the  assumption  of  the  independence  of  particles  1  and  2 
from  the  rest  of  the  system  in  the  lowest  order  approximation  becomes  in¬ 
valid  if  the  particles  have  an  infinite  range  of  correlation.  The  presence 
of  a  third  particle  cannot  be  ignored. 

If  there  is  no  collision  between  the  time  0  and  t  ,  S 

t  %  /  .  \  O 

and  coincide  from  x  *  0  and  on  to  infinite,  the  integral 

is  zero  identically* 


Next  we  consider  the  contribution  of  the  integral  terms  to  the 
solution  for  g^(x, 12,0 .  Since  the  structure  of  the  second  square 
bracket  of  the  integral  term  in  Eq.  (12)  is  very  much  the  same  as  the 
first  one,  we  shall  analyze  the  first  one  only,  i.e., 


n 

m 


d—  d— 3 


d$ 

d£ 


37*  f s_t  3,t)  -  S_a>(C,13)]f(o)(l)f(o,(?)f(o,(3) 


(15>) 


It  is  seen  that  because  of  its  being  x  and  t  indeDendent.  the 

—  o  •  ’ 

second  term  in  the  above  expression  will  cause  a  secularity  in  the  solution 
for  To  be  free  of  secularity  one  usually,  though  not  alwavs  cor¬ 

rectly,  seeks  a  cancellation  between  the  first  and  second  terms  for  all 
£  and  £. .  This  led  Sandri^  to  conclude  that  the  cancellation  is  not 

exact  for  some  kinds  of  interaction  between  the  three  particles  and  thus 

£ 

he  was  faced  with  a  general  secularity  (non-local)  in  the  exnansion. 

However,  considering  the  expression  in  Kq.  ( 1 5 )  as  a  source  function  of 

the  differential  equation  for  g^-(x,12,t  ),  we  see  that  the  solution 

for  g^^(x_,l?,to  •*  »)  depends  very  much  on  the  behavior  of  this  source 

function  in  terms  of  x.  It  is  evident  that  if  this  source  function  has 

a  finite  range  in  £,  and  goes  to  zero  sufficiently  fast  as  |xj  ■*  °°,  the 

solution  for  g^^(x,12,t  )  will  be  well-behaved  in  the  limit  of  t 
~  o  o 

To  investigate  the  variation  of  the  source  function  in  terms  of 
3t,  we  consider  that  at  time  t  the  configuration  of  particles  1,  2, 
and  3  is  as  shown  in  Fig.  2.  We  assume  that  the  three  particles  inter¬ 
act  simultaneously  at  the  time  X  .  Apart  from  the  integral  operator 


Sandri  used  the  hierarchy  equations  for  S-particle  functions;  the 
corresponding  quantity  which  caused  the  secularity  is  then 


•r3 


where  1',  and  3'  are  the  velocities  of  the  three  particles  at  the 

ooints  -t  backwards  along  the  particle  trajectories  from  the  initial 
o 

configuration  123.  For  the  second  term  in  the  source  function  the  oarticle 
;>  does  not  interact  with  either  1  or  3.  The  trajectories  of  1  and  3  are  ’ 
then  given  by  the  dashed  lines  in  Fig.  2  because  of  the  absence  of  2. 

Thus 

f(o)(2)  f.  (4,13)  =  f(o)(?)  f(o)(l")  f(o)(3")  (17) 

•tOO  —  ’ 

The  points  l"  and  3"  should  be  projected  infinitely  far  back  from  1  and  3. 
After  particles  1  and  3  leave  the  interaction  zone,  there  is  no  change  in 
their  velocity;  it  is  immaterial  exactly  where  l"  and  3"  are,  orovided 
they  are  out  of  the  interaction  zone.  For  such  an  event,  (16)  and  (17) 
are  not  the  same  in  the  limit  of  t^  -*■  however  their  difference  is 
finite.  As  |x|  increases,  the  effect  of  the  particle  2  on  the  motion 
of  particles  1  and  3  becomes  less.  For  a  sufficiently  large  |xj  we  see 
that 

1'  -  1" 

2'  ♦  2" 

3’  -  3" 


and  the  expressions  in  (16)  and  (17)  exactly  cancel  each  other  except 


in  the  following  two  instances 


action  1*2*3.  The  source  function  does  not  vanish.  It  can  be- seen  i 


that  the  expressions  (l6)  and  £l?)  do  not  cancel  each  other  even  in  the  ! 

limit  of  |xj  •*  »  (see  Fig.  3).  However  in  such  a  limit  we  see  that 
I -  *|  also  tends  to  be  very  large.  This  follows  because  |sj  has  to 

be  finite  because  of  the  factor  dc/d£_  in  the  source  function.  It  can  be  .4 

-  /i 

.4 

seen  that  in  the  limit  of  - _x|  ■+  "  in  order  that  the  second  binary 

collision  3*2  occur  at  all,  the  relative  velocity  between  2  and  3  must  be  ;i 

selected  very  carefully,  'lhe  relative  velocity  between  3  and  2  after  ‘5 

■J 

1*3  (oarticles  3  at  3')  must  be  oriented  almost  in  the  direction  of  ?'i 

•  o 

-  x.  In  other  words,  particle  3  at  3^  must  aim  at  particle  2  which 
is  at  a  large  distance  away.  The  solid  angle  in  the  relative  velocity 


space  ,  within  which  the  second  binary  collision  is  possible,  goes 
o 

as  1/|C,  -  *1  ?  as  |£  -  x|  -►  ”,  We  thus  conclude  that  although  the  ex¬ 
pressions  in  (16)  and  (17)  do  not  cancel  in  the  limit  as  |x|  -*■  ”,  the 
source  function  itself,  which  contains  a  velocity  integration  dv^+,  Poes 
to  zero  in  the  limit  as  |xj 

2)  Two  successive  binary  collisions  1*3  and  1*2.  The  situation 
is  much  the  same  as  before.  In  order  to  have  such  successive  binary 


collisions  in  the  limit  as  |x|  ^  we  require  a  special  orientation 
of  the  relative  velocity  v^,  Since  now  there  is  no  integration  over 
either  or  Vg,  the  source  function  for  this  special  arrangement  (or 
equivalently  the  special  local  region  in  the  phase  space)  becomes  in- 


The  symbol  *  between  two  numbers  indicates  the  interaction  between 
the  two  particles  designated  by  the  corresponding  number. 

+  Note  velocity  v^,  and  v^  are  related  in  a  definite  way.  The  integration 
in  v*  amounts  to  an°integration  in  v_,, 

mm~-  ’•  -i 1  -A  • 


■  :1 


finite  in  range.  The  contribution  to  the  solution  for  g^(x,12,to -*■  “ ) 
becomes  secular.  This  local  secularitv  is  very  similar  to  that  caused 
by  the  term  dg^°\/dct.  However  these  two  singular  terms  do  not  in 
general  cancel  each  other  exactly.  In  order  to  ensure  the  secularity 
of  Eq.  (13),  the  orientation  of  v^  is  determined  indenendently  of  the 
presence  of  narticle  3,  while  in  Eq.  ( 15 )  there  is  a  strong  dependence 
on  narticle  3  (comnare  Fig.  lb  and  Fig.  3).  We  therefore  conclude  that 
the  solution  for  g^ 1  ^ (3c,l?,t )  is  locally  secular  at 


x  =(—  (v  -►  3)  t.  -► 

—  i  —it  —  * 


(16) 


This  exDression  includes  Eq.  (lM  as  a  srecial  case. 

We  have  indicated  earlier  that  in  this  soecial  local  region, 
g^°^(jc,12)  is  non-secular  but  has  an  infinite  range  in  x,  Now 
g^^(x,12)  becomes  proportional  to  tQ  as  tQ  becomes  large.  In  order 
to  preserve  the  asymptotic  character  of  our  series  representation  of 
g  in  the  limit  as  t  -*  we  must  reouire 

t«i 

e 


(19) 


and  consequently  from  Eq.  (18)  and  (19) 


lil  «  T 


(20) 


Since  our  time  and  length  are  normalized  by  the  collision  time  and 
range  of  the  pair  potential  respectively,  the  inequalities  (19)  and  (20) 
indicate  that  the  formal  expansion  so  far  carried  out  is  at  most  valid 
only  for  a  time  duration  of  the  order  of  the  mean  free  time  and  a  length 
scale  of  the  order  of  the  mean  free  path.  The  solution  for  g^  cannot 
be  valid  in  this  special  region,  since  the  binary  interaction  alone 


cannot  be  the  leading  approximation.  In  other  words,  when  we  solve  for 
s(x,|12)  (the  correlation  function  of  1  and  2),  the  presence  of  a  third 
particle  can  no  longer  be  ignored. 


We  now  introduce  the  following  scaling  for  the  two-particle  cor¬ 


relation  function  equation: 


(21) 


'J 


where  t'  and  x-  will  be  of  order  unity  for  a  time  of  the  order  of  the 
mean  free  time  and  a  particle  separation  of  the  order  of  the  mean  free 
path.  Using  Eq.  (21)  in  Eq.  (2),  we  obtain 

at-  +^12  *  If-  -  iff.  *  [e(12)  +  f{1)f(2>i  <22> 

—  —  — 1  —2 

* 

*  “  d£  dv^  •  (j^-  [  f(l)gU-x;23)  ♦  h(x'  ,c.,123)  ] 

-  (f(2)  g(x»-£,13)  «•  h(x’,x'-£,123)]) 

This  equation  cannot  be  solved  independently  of  f  as  before,  since  both 
f  and  g  are  now  t-  dependent.  We  see  that  the  governing  system  in  this 
local  region  is  non-Markoff ian  and  Bogoliubov’s  functional  assumption 
is  not  valid.  In  other  words,  the  evolution  of  the  system  in  this  local 
region  cannot  be  described  by  a  self-contained  one-particle  equation; 
the  correlation  functions  vary  in  the  slow  time  scale  on  their  own  right 
(not  as  functional  of  the  one-particle  function). 


Kovever,  the  contribution  of  such  a  local  nonMarkoffian  solution 
to  the  one-particle  function  is  only  a  higher  order  effect  because  of 


For  a  pair  potential  having  a  range  much  smaller  than  the  mean  free 
path,  the  contribution  of  g  to  f  from  the  above  mentioned  local  region 
is  negligibly  small.  We  therefore  conclude  that  as  far  as  the  one- 
particle  function  is  concerned,  one  can  ignore  completely  the  local 
non-Markoffian  region  indicated  earlier.  In  other  words,  for  the 
correction  term  of  the  Boltzmann  collision  integral,  one  simple  solves 
Eq.  (12)  for  g^^t  •*  «)  and  substitutes  it  into  En,  (23). 

The  solution  for  g^(x,12,t  )  is  obtained  from  Eq.  (12)  as  follows: 

-  o 

g(l)(x,12,to)  +  f(l)(l,tQ)  f(o)(2)  +  f(1)(2,to)  f(o)(l)  (2U) 


the  velocity  integration  dv^  and  the  potential  factor  d$/dx  in  the 
one-particle  equation,  i.e.. 


dx  dv  ^ 
—  —2  dx 


*1 


(23) 


=  C_t  (12)  [f(l)(l,to=0,t1)  f(°)(?,t1) 

O 


+  f(l)(2,to=0,t1)  f(o)(l,t1) ] 


♦  a 

m 


dT  S  (12) 

-T 


— 1  o 


-S 


-  =-t  .,(12>  IT"  -  ITT  s-t 

o  — 1  —2  o 


♦  5  ^T(12)  jf-  S_J-£,23))  f(o)(l)f(o)(2)f(o)(3) 

*  o  d— 2 


where  f^Vt  »0)  is  the  initial  value  of  f^(t  ). 

o  — # 

For  the  investigation  of  the  one-particle  equation,  the  second  and 


is 


j-  -til  'tie*  *3$,*  * £ 


1  r*»  i  !►*  '  *  /  _  ^  -  --  !» 


third  terms  on  the  left  side  of  Eq •  (2k)  can  be  ignored  because  of  the 

3C  integration  in  Eq,  (23).  Given  a  finite  x,  for  sufficiently  large,  t  , 

the  variation  in  t  of  the  first  two  terms  inside  the  curly  brackets  of 

Eq.  (2k)  is  as  shown  in  Fig.  k.  The  little  bump  hear  t  *  t  represents 

o 

the  interaction  of  the  particles  at  the  time  t  ■  Mote  that  the  successive 

o 

binary  collisions  discussed  earlier  have  been  neglected.  In  other  words, 
the  local  region  as  indicated  in  Eq,  (l8)  is  excluded  in  the  present 
discussion. 

The  operator  S_^(l2)  makes  the  contribution  of  a  function  which 
has  a  finite  range  in  |x|  ,  completely  negligible  in  the  limit  of  large 
Therefore  as  far  as  the  integral  over  t  in  Eq.  (2k)  is  concerned,  we  can 
neglect  the  bump  in  Fig.  k  and  the  corresponding  expression  for  Eq.  (2k)  is 

F^  (x,12,to-  -,to)  =  g(l)(l2)+f(1,(l)f(o)(2)4.f(l)(2)f(o)(l) 


S_^(12)  [f(l)(l,to)f(o)(2,t1)  ♦  f(l)(?,t1)f(o)(l,t1)J 


+  |  dT  S_x<12>  di  di3  *  {3“  S^.(x,i.l?3) 


S-U2)  S-(i’13)  ‘  ^ 


-  S_t<(12)  |~S_Ji,23))  f(o)U)f(o)(2)f{o)(3) 


This  expression  for  Fg  '  is  equivalent  to  the  one  obtained  by  Choh 
and  Uhlenbeck  except  for  the  first  term  on  the  right  side,  which  vanishes 
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for  f^(t  =0)  =  0.  It  was  first  uointed  out  by  Sandri^  that  the  role 
of  the  functional  dependence  in  the  Bogoliubov  method  was  played  by 
the  term  3g^°^/3t^  in  the  multiple-time-scale  formulation.  Still  ignoring 
the  successive  binary  collisions  discussed  earlier,  the  solution  of  Eq. 
(25)  can  be  written  in  the  following  more  familiar  form: 

F  (x,12,t  -  »)  =  S  (12)  [f(l)(l)f(o)(2)  +  f(l^(2)f(o)(l)] 
d  "*  o  — *» 


d£  dv.  [S  (x, 5,123)  -  S  (x,12 )  S  (5,13) 

*  “j  _QO  — “  —  aiOD  aOO 


(26) 


*  S^tx,!?)  S_oj(5-x, 23)  +  S_<_(x,12)  ]  f(o)(l)f(o)(2)f(o)(3) 


The  corresponding  one-particle  equation,  ignoring  the  transient 
of  order  of  the  collision  time,  is  as  follows: 


d£dv  S  (x.^Hf^M^+f^M0^] 

- 3  -°°  — 

di  di3  (S_J*.£»123>  (27) 

-  S_-(x,12)  5_j£,13)  -  S_Jx,12)  S_JS.-x.23) 

♦  Sjx.12)]  f(o)(l)f(o)(2)f(o)(3) 

The  variation  of  f^  on  the  t,,  time  scale  can  only  be  determined 
by  the  condition  of  non-secularity  of  f^  in  the  long  time  limit  of  t^. 

In  the  limit  of  t^  -*  »,  by  the  well  known  H-theorm  of  the  Boltzmann 
(binary)  collision  integral,  f^  goes  to  the  Maxwellian  distribution. 


n 

m 


dx  dv 


-  -2  dx 


.  JL 


i 


3f 


(1) 


at. 


3f(o)  n 
3t_  ”  m 


dxdv  ^ — 

- 2  ^ 


It  is  not  difficult  to  show  that  the  right  hand  side  of  Eq.  (27)  vanishes 


equal  to  the  Maxwellian  distribution.  Therefore  we  have 


the  condition  of  non-secularity ,  that 


(28) 


3fll)  n  [  .  .  dt  3 

— rr -  ■  —  dx  qv  — *•  •  - — - 

3t.  b  '  —  — 2  dx  3v, 

1  j  -  -1 


fv  *  -dr  s  Ji2)  (f(1)(i)f(o)(2)+f(1)(2)f(o)(i) 


♦  |  di2  «  •  3^  di  %  lS-(i*i'123)  - 


-  S_Jl2)  S_oo(£-x,23)  ♦  S_a>(x*12)]f(o)(l)f(o)(2)f(o)(3)  (29) 


The  second  terra  on  the  right  side  of  Eq.  (29)  is  the  triple  collision 

integral  which  has  been  obtained  by  Choh  and  Uhlenbeck.,^  Oreen,^  Rice, 

7  8 

Kirkwood,  and  Harris,  and  Resibois.  Equation  (20)  indicates  a  self- 
termination  of  the  multiple-time-scale  formulation,  i.e. ,  there  is  no 
variation  after  the  time  scale  of  kinetic  evolution,  which  is  what  one 
would  expect  physically  for  a  homogeneous  system. 

It  is  well  known  that  on  the  binary  collision  level  a  system  with  any 
disturbance  will  go  to  the  thermal  equilibrium  as  indicated  by  the  H- 
theorem  of  Boltzmann.  It  is  of  interest  to  see  whether  Eq.  (29),  to¬ 
gether  with  Eq.  (9),  will  force  a  system  into  thermal  equilibrium,  i.e., 
the  one-particle  function  becomes  the  Maxwell-Bolt zmann  distribution  up 
to  the  order  of  c.  Since  the  expansion  of  the  hierarchy  in  terms  of  c 
is  viewed  as  an  asymptotic  series,  it  is  only  consistent  to  answer  the 


above  nuestion  order  by  order.  First,  because' of  •  the  H-theorem  of  the 
dinar-'  collision  integral  we  know  becomes  Vaxwelliar.  in  thenr.nl 

(2] 

equilibrium.  '.lith  f  =  f  ,  the  triple  collision  integral  vanishes 

max 


and  we  are  left  with 


iif 


(1) 


it. 


dxdv. 


iil  .  _L_  c- 

dx_  3-v  - 


(1?) 


f( 1 ) ( 1 ) f < ( ? )+f C 1 )  ( p ) f C 0 ) ( 1 )  ] 


(jol 


Combining  this  with  the  lowest  order  approximation  (Boltzmann's  equation), 
we  obtain 


-f-(f(o)  ♦  Ef(l))=^ 

m 


dxdv 
- ?  dx 


aii 


f(o)(D+cr(1)(i): 


]  t  r(  0  J  ( p  1  ^  ( 2 )  )+o(e2) 


By  the  H-theoreir.  for  this  collision  integral,  we  know  that  f^°'  +  ef^^ 
goes  to  the  Maxwellian  distribution  in  thermal  equilibrium.  If  we  choose 
f^)(to=0)  =  0  initially,  f^^  poes  to  zero  in  thermnl  equilibrium  because 
Eos.  (P 9)  and  (30)  preserve  the  normalization  of  f'  .  i'uch  a  statement 
apparently  holds  also  for  the  higher  order  approximations. 


IV.  Discussion. 

Using  the  multiple-time-scale  method,  it  is  found  that  the  triple 
collision  integral  for  a  classical  3oltzmann  gas  is  identical  to  that 
Obtained  by  Choh  and  Uhlenbeck,  The  functional  deoendence  of  the 
multiple  particle  functions  on  the  one  narticle  function  is  nevertheless 
found  to  be  locally  invalid  for  particles  having  large  separations  at 
time  tQ  and  with  their  relative  velocity  oriented  in  such  a  way  that  they 
were  in  collision  at  t^  *  0.  Such  locally  non- functional  behavior  of 
the  multiple  particle  functions  is  related  to  the  inseparability  of  the 
time  .and  length  scales  caused  by  the  two  successive  binary  collisions 
which  happen  at  a  time  interval  of  the  order  of  the  mean  free  time.  A  new 
scaling  for  such  a  local  region  is  given.  The  self-contained  (or  Markoffian) 
one-particle  description  of  the  system  which  is  true  outside  such  a  region 
is  found  to  be  impossible.  The  hierarchy  equations  have  to  be  solved 
simultaneously  on  the  kinetic  time  scale  (of  the  order  of  the  mean  free 
time).  However,  due  to  the  very  special  nature  of  the  relative  velocity 
in  this  local  breakdown,  the  contribution  of  this  local  region  to  the 
one-particle  function  comes  into  play  only  at  orders  higher  than  the 

triple  collision  level.  _ Furthermore ,  if  the  pair  potential  is  limited 

* 

in  its  range,  3uch  a  local  contribution  can  be  neglected  even  in  higher 
orders. 
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